Shear viscosity to entropy density ratio in nuclear multifragmentation 
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Nuclear multifragmentation in intermediate energy heavy ion collisions has long been associated 
with liquid-gas phase transition. We calculate the shear viscosity to entropy density ratio rj/s for an 
equilibrated system of nucleons and fragments produced in multifragmentation within an extended 
statistical multifragmentation model. The temperature dependence of rj/s exhibits surprisingly 
similar behavior as that for H2O. In the coexistence phase of fragments and light particles, the 
ratio -q/s reaches a minimum of comparable depth as that for water in the vicinity of the critical 
temperature for liquid-gas phase transition. The effects of freeze-out volume and surface symmetry 
energy on rj/s in multifragmentation are studied. 
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Understanding the behavior of nuclear matter under 
extreme conditions has been one of the most important 
present challenges in heavy ion physics. Multifragmenta- 
tion in intermediate energy heavy ion collisions provide 
a key mechanism to address this issue where an excited 
nucleus formed expands and breaks up into various frag- 
ments and light particles [1H3- The final yield distribu- 
tion is quite sensitive to the internal excitation, breakup 
density of the nucleus and the symmetry energy part of 
the binding energy of the fragments 0-13 • 

Due to van der Waals nature of the nucleon-nucleon 
interaction, it is expected that multifragmentation may 
exhibit features of liquid-gas phase transition [l|, |2| . Evi- 
dence of this provided from the observation of the nuclear 
caloric curve relating excitation energy to the tempera- 
ture of the breakup source |9l4ll[. Extensive studies have 
been carried out to understand the dependence of the 
caloric curve on the system size and particularly on the 
density dependence of nuclear symmetry energy which is 
poorly constrained jl2| . 

Multifragmentation studies are however mostly con- 
fined to the effects of the state variables on the thermo- 
dynamic properties of the system. Whereas transport 
properties of the dynamically evolving system of nucle- 
ons and fragments formed in fragmentation have received 
little attention. Since the transport coefficients char- 
acterize the dynamics of fluctuation of the dissipative 
fluxes in a medium 13|, their knowledge is essential for 
a better understanding of the fragmentation observables. 
Recently, considerable attention has been focussed pri- 
marily on the shear viscosity coefficient that involves the 
transport of momentum due to velocity gradient in an 
anisotropic medium. Empirical observation J14l | of the 
temperature dependence of the shear viscosity to entropy 
density ratio rj/s for water (as well as He, Ne2) exhibits a 
minimum in the vicinity of the critical temperature Tc for 
liquid-gas phase transition. Furthermore, a lower bound 
of r]/s > l/in obtained by Kovtun-Son-Starinets (KSS) 
[l5{ in certain gauge theories is speculated 15[ to be valid 
for several substances in nature. 



In heavy ion collisions at energy E^iab ^ 1000 
MeV/nucleon, the shear viscosity coefficient has been es- 
timated for nucleon transport in the Uhlenbeck-Uehling 
equations flq. Analysis of the observed transverse flow 
of nucleons in the microscopic QMD model in Au-|-Au, 
Nb-fNb and Ca+Ca central collisions at energy _Eiab — 
400 — 1200 MeV/nucleons requires shear viscosity of 



30 < ?7 < 60 MeV/(fm'^c) 



Since the degrees of 



freedom in these studies involve only nucleons hence the 
r] estimate is relevant only in the early stages of the reac- 
tion. Whereas fragments formed due to correlations and 
fluctuations of the system that goes beyond the mean 
field dynamics of nucleon transport should result in a 
different shear viscosity. At freeze-out the momentum 
transport and thereby the shear viscosity is mostly ef- 
fected by long range Coulomb interaction between the 
charged particles. 

In this paper we estimate the shear viscosity to en- 
tropy density ratio rj/s of an equilibrated system of frag- 
ments and/or light particles produced from multifrag- 
mentation within a modified microcanonical statistical 
multifragmentation model (SMM) [21 that has success- 
fully reproduced several observables |2h4|. I18|. We shall 
show that rj/s decreases as a function of temperature of 
the system, exhibits a minimum at the coexisting liquid- 
gas phase comprising of intermediate mass fragments and 
light particles, and increases again at high temperature 
in a system of light particles. The estimated minimum of 
rj/s is found comparable to that for water in the vicinity 
of critical point for liquid-gas phase transition. 

In the SMM a hot source with mass and charge 
{Aq, Zq) at temperature T expands to a freeze-out volume 
y = (l + x)l^ and undergoes prompt statistical breakup; 
here Vo is the normal nuclear volume and x > is input 
parameter. We start with the grand-canonical version of 
SMM which consists of minimizing the free energy F of 
the system: 
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under conservation of mass Aq — "^a z ^A,zA and 
charge Zq — ^^ ^ N^zZ- The translational energy of a 
nuclear species with mass and charge (A, Z) is 
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where the thermal wavelength is At = y^nh /mT, m 
is the nucleon mass. The spin degeneracy is denoted by 
gA,z and the free volume Vf =V — Vq. The multiplicity 
Na,z of the fragment {A, Z) is then given by [2, |8[ 
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The baryon and charge chemical potentials /is and [iq 
are obtained from mass and charge conservation. The 
internal free energy of the species (A, Z) is 
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where the parameters are eo — 16 MeV, /3o — 18 MeV 
and Tc = 18 MeV. The fragments here are assumed to 
be at normal density. The Coulomb repulsion between 
the fragment F^ in Eq. ([T]) is evaluated in the Wigner- 
Seitz approximation and corresponds to the last term in 
Eq. (H]). The Coulomb self energy of a fragment is in- 
cluded in its binding energy Ba,z- The last two terms in 
Eq. ([3]) stem from homogeneous term of the Wigner-Seitz 
approximation [2|, |8| . Light nuclei with A <b are consid- 
ered point particles with no internal degrees of freedom 
so that the bulk and surface contributions to internal free 
energy /^ ^ [second and third term of Eq. (|3]), respec- 
tively] are neglected, except for alpha particle, where the 
bulk contribution is retained as usual [^ . For these light 
nuclei, experimental values for binding energy are used. 
For heavier nuclei A > 5, the spin degeneracy gA,z = li 
and for binding energy Ba.z, the computed Liquid Drop 
Mass (LDM) formula of Ref. is used: 



Ba,z = CvA - CA^I^ -C,^^ + Cd 
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where Ci = ai\V — ki{A — 2ZY/A^] and i — v,s cor- 
responds to the volume and surface terms, respectively. 
The volume and surface contributions to symmetry en- 
ergy give E'sym = Csym{A - 2Z)^/A'^ where Csym = 
ayky—asks/A^/^. Two versions of LDM formula are used 
here, dubbed as LDMl and LDM2 [8], to study surface 
symmetry energy effects on the rj/s ratio. The simpler 
LDMl has fc^ = (also Cd = 0) and thus neglects the 



surface corrections to the symmetry energy. The com- 
plete LDM2 formula preserves all the terms. 

The energy conservation of the fragmenting source 
leads to 

Ef^^^ + E*= E'^{T) + ^ Na,z [-Ba,z + eXz] 
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where Ef^^^ is the ground state energy of the source 
and Ea z (^) denotes the excitation energy of the frag- 
ment {A, Z) at temperature T with a translational energy 
Ea z — 3r/2. Equation © allows to extract the excita- 
tion energy E* of the source at a given T and thereby 
construct the nuclear caloric curve. 

In the grand-canonical SMM the probability distribu- 
tion of the fragment yield, Pa,z = Na,z/ X^a z ^a,z, al- 
lows one to generate [19[ a Monte Carlo microcanonical 
ensemble of fragments with exact conservation of mass 
Aq, charge Zq and total energy of the source. The frag- 
ments in a microcanonical ensemble are placed in a non- 
overlapping fashion within a spherical freeze-out volume 
V. These particle are then allowed to evolve in time 
under Coulomb repulsion within the freeze-out volume 
with periodic boundary conditions in the configuration 
space. Finally, collision between the fragments enforce ki- 
netic equilibration when the system is found to approach 
momentum isotropization [20|. The shear viscosity of 
this dynamically evolving equilibrated system is then es- 
timated. 

In Fig. [T] we show the charge distribution for the 
breakup of i^OSm at T = 3,4,6 MeV in the LDMl and 
LDM2 mass formulas. At low temperature T = 3 MeV 
the system is characterized by few light particles and a 
massive nucleus (liquid phase) close to the source charge. 
The minimum of the free energy F = E — TS is es- 
sentially controlled by the surface term in Eq. @ that 
favors one massive nucleus instead of small fragments 
with a larger total surface. With increasing temperature 
more intermediate mass fragments, IMFs {3 < Z < 15), 
are produced. Here the —TS term in free energy, domi- 
nated by binding energy and Coulomb repulsion, favors 
the breakup into small sized fragments [3l. Eventually 
at large T only light particles are produced (gas phase). 
Compared to the LDMl set, the inclusion of surface sym- 
metry energy in LDM2 that are important for light nuclei 
of mass A < 25 |8i|, leads to the suppression of isospin 
symmetric IMFs and an enhancement of neutron rich 
heavier nuclei. 

In order to get insight into the quantitative effects from 
the two mass formulas, we show in Fig. [21 the multiplicity 
of different species as a function of temperature. At low 
T < 3 MeV where the yield is dominated by a massive 
fragment, the symmetry energy effects are imperceptible 
in the LDMl (triangles) and LDM2 (solid circles) sets at 




FIG. 1: (Color online) Charge distribution in the breakup 
of ^'' Sm nucleus at various temperatures T at a freeze-out 
volume oi V = 6Vo in the LDMl (thin hues) and LDM2 
(thick lines) mass formulas. 



the same freeze-out volume V — GVq. At moderate and 
high temperatures the total multiplicity iVtot in LDM2 
set is reduced as surface symmetry energy favors heavier 
nuclei with large neutron-proton asymmetry. This leads 
to a reduction in the yields of more isospin symmetric 
IMFs, iViMF, and light nuclei, A^ip, in comparison to the 
LDMl set. Also shown in Fig. [2]the particle abundances 
in the LDM2 set at a smaller freeze-out volume V = 3Vb, 
where as expected, the yields are reduced. 
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FIG. 2: (Color online) The temperature dependence of mul- 
tiplicity of (a) all particles, A'tot, (b) IMFs, A^imf, (c) light 
particles with A < 5 except alpha, Nip, and (d) alpha, Na, in 
the breakup of ^''"Sm nucleus in LDMl (triangles) and LDM2 
(solid circles) mass formula at a freeze-out volume oiV — 6Vb 
and in LDM2 (open circles) at 1/ = 3Vb. 



The nuclear caloric curve |2l.l9l4ll associated with the 
fragmentation of ^^"^Sm nucleus is shown in Fig. [31 At 
3 ^ £'*/A < 8 MeV, the caloric curve exhibits a slow 
but monotonous increase of temperature with excitation 
energy [10|, |2l|, \2^ which may be a signature of liquid-gas 
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FIG. 3: (Color online) Caloric curve expressing excitation en- 
ergy per nucleon E* /A versus temperature T for the breakup 
of ^^°Sm nucleus in LDMl (triangles) and LDM2 (solid cir- 
cles) mass formula at a freeze-out volume of V = 6Vb and in 
LDM2 (open circles) at V' = SVq. 



phase transition. This behavior stems from energy con- 
servation constraint of Eq. ^ when appreciable amount 
of energy is used to produce abundant IMFs and light 
particles. While the experimentally observed plateau 
[lOl lllj l in the caloric curve may result if the breakup 
occurs at a fixed pressure (Le. multifragmentation is 
an isobaric process) [2l|, [23 



24| in contrast to the fixed 



freeze-out volume employed here, or if the fragments are 
expanded 25| unlike the fragments assumed here to be at 
normal nuclear density. Note in this excitation range, the 
massive fragments with large binding and internal excita- 
tion in the LDM2 set (solid circles) lead to slightly higher 
temperatures (or conversely smaller E* / A) compared to 
LDMl set (triangles). While at high E*/A> ^ MeV, 
the smaller multiplicity of light particles, which have no 
internal degrees of freedom, in the LDM2 set results in 
lower temperatures relative to LDMl set. At freeze-out 
volume V = 3Vb in LDM2, the breakup temperature is 
consistently higher as more massive and fewer fragments 
are produced. 

The entropy S = —dF/dT in the microcanonical en- 
semble of fragments is determined using the conventional 
thermodynamic relation [2, |3| 
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where the total multiplicity in an event is M ^ 
^A z ^A,z and the last term corresponds to entropy con- 
tribution from the bulk and surface terms in the internal 
free energy of Eq. ^ . 

In Fig. mja) the entropy per nucleon S/A is shown as 
a function of temperature for fragmentation of ^'"^"Sm. In 
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FIG. 4: (Color online) (a) Entropy per nucleon S/A, and (b) 
shear viscosity coefficient 77 as a function of temperature for 
the breakup of ^^"Sm nucleus in LDMl (triangles) and LDM2 
(solid circles) mass formula at a freeze-out volume of V' = 6Vb 
and in LDM2 (open circles) at 1/ = 3Vo. Shear viscosity in 
nucleon transport calculation [Ig is shown in diamond. 



the temperature range 4 < T < 6 MeV associated with 
liquid-gas mixed phase in the caloric curve, the entropy 
shows a rapid increase with temperature. In this region 
abundant IMFs are produced in the LDMl set (trian- 
gles) with smaller internal degrees of freedom that leads 
to somewhat higher entropy than in LDM2 set (solid 
circles) with more massive fragments. At I^ = 3Vo in 
LDM2, suppression of particle yield leads to a reduction 
in entropy. 

For time evolving system of fragments in equilibrium, 
the shear viscosity due to momentum transport (via 
Coulomb scattering between the fragments) can be com- 
puted from Kubo relation [20, '2&\ or in the classical ki- 
netic theory [27[. The Kubo formula employs the lin- 
ear response theory to relate the transport coefficients 
as correlations of dissipative fluxes. However it provides 
the total viscosity of the system and not from individ- 
ual species. On the other hand, in the kinetic theory, 
the total shear viscosity of a multicomponent system can 
be expressed as the sum from individual contribution as 
rj c± (1/3) X^i '^j(pj)-^i; where Ui is the number density, 
(pi) is the average momentum and A,; is the mean free 
path of the zth species. Moreover, A^ = 1/ J2j '^j^ij 
where cr,,- is the coUisional cross section which is taken as 



that for the usual Coulomb scattering. The average ther- 
mal momentum of the particle in the nonrelativistic limit 
is {pi) — mi{vi) — ^J^mi{rlTx. The results presented be- 
low are in the kinetic theory limit; we have checked the 
total T] so obtained matches with that from the Kubo 
formalism. 

In Fig. m^b) we present the shear viscosity 77 for an 
equilibrated ensemble from fragmentation of ^^'^Sm. At 
temperatures T < 3 MeV (corresponding to the liquid 
phase), the shear viscosity is found to rapidly increase 



with decreasing T. For the dilute system comprising of a 
large nucleus and few light particles the mean free path 
A is large. Thus a fragment can transport momentum 
over a large distance resulting in large r\. At interme- 
diate T ^ 3 — 6 MeV (corresponding to the coexistence 
region) , the magnitude of viscosity is determined by com- 
peting effects between the size and multiplicity of the 
fragments in the system. The Coulomb repulsion will 
force the two colliding nuclei to occupy the available free 
space (void) which in turn will collide with the neigh- 
boring nucleus and so on. This procedure can effectively 
transport momentum over a large distance and produce 
a large viscosity [IJ]; the coefficient grows with temper- 
ature as 77 '^ T^l"^ . In general, the viscosity r\A,z of a 
species (^, Z) was found to progressively increase from 
heavier to lighter particles that have larger ZjA ratios. 
Compared to LDMl ?XV — 6Vb, the heavier fragments 
in LDM2 in the liquid-gas mixed phase and somewhat 
higher temperature (see also Figs. [2]and[3l) are more ef- 
fective for momentum transport resulting in an increased 
77. This is due to r\A,z values are ~ 30% higher in the 
LDM2 than in LDMl at a given T . At a smaller freeze- 
out volume Y = 3Vo fewer fragment collisions as well 
as smaller density of voids inhibits momentum transport 
and thereby reduces dissipation in the medium. 

For orientation we also show in Fig. IH^b) the shear vis- 
cosity from analytical fit to the numerical evaluation of 77, 
obtained by Danielewicz [16], in the Uhlenbeck-Uehling 
transport equations within a first order Chapman-Enskog 
approximation: 



77 = (1700/r2)(n/no)2 + [22/(1 

+5.8ri/V(i + leoT^) . 



(8) 



In this microscopic calculation, where the relevant de- 
grees of freedom are nucleons in the colliding nucleus, sig- 
nificant momentum transport by the nucleons in the early 
stages of reactions (not considered in our study) leads to 
faster growth in viscosity at small T ^ 1 — 2 MeV. In 
contrast, at T > 4 MeV, the present study clearly under- 
scores the importance of finite size fragments from multi- 
fragmentation (missing in the transport calculations [16]) 
that substantially enhance the viscosity in the medium. 
Figure [5] shows the shear viscosity to entropy density 
ratio 77/s as a function of temperature in multifragmen- 
tation of ^^"Sm. The 77/s value gradually decreases as 
a function of rising temperature up to T ~ 4 MeV in 
the liquid phase (system dominated by a massive nu- 
cleus) and then increases again at T > 6 MeV in the gas 
phase (system of light particles) . At T « 5 MeV that is 
close to the critical temperature for multifragmentation 
(liquid-gas coexistence phase) [l|-y, Isj the 77/s reaches a 
minimum. In the LDM2 mass formula at a freeze-out 
V = 6V01 a- minimum of (77/s)niin ~ 2.1 is obtained. 
In the absence of surface symmetry energy in LDMl, 
the (77/s)inin turns to be somewhat smaller. While at a 
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FIG. 5: (Color online) The shear viscosity to entropy density 
ratio r]/s as a function of temperature T for the breakup of 
^^°Sm nucleus in LDMl (triangles) and LDM2 (solid circles) 
mass formula at a freeze-out volume of 1/ = 6Vo and in LDM2 
(open circles) ai V = 3Vo- The rj/s from nucleon transport is 
shown in diamond. The inset shows ri/s vs T for water from 
Ref. [Ijl for an isobar at the critical pressure, Pc = 22.6 MPa, 
(dashed lines) and one below it at P = 10 MPa (solid line) 
and the other above it at P = 100 MPa (dotted line). 



smaller freeze-out volume oi V — 3Vb, the ij/s magni- 
tude is found to be even smaller. The (?7/s)min estimated 
in multifragmentation is however significantly above the 
conjectured KSS lower bound of l/in 15,]. If we adopt 
the entropy of the LDM2 set with V — GVq and rj of Eq. 
dS]) for nucleon transport [1^1, then 77/s is seen (Fig. [5]) 
to decrease continuously with increasing T. 

Interestingly, the (??/s)min obtained in multifragmen- 
tation is comparable to the minimum value for isobars 
passing in the vicinity of the critical temperature Tc for 
liquid-gas phase transition in water [IJ]. For H2O when 
an isobar passes through the critical point (shown in the 
inset of Fig. [5]), the (??/s)min forms a cusp at Tc. When 
the isobar passes below Pc the minimum is at T < Tc with 
a discontinuous change across the phase transition. For 
an isobar passing above Pc, a broad minimum is found 
at a T slightly above Tc- In fact the smallest value of 
r]/s corresponds to the most difficult condition for the 
transport of momentum. In analogy to this observation, 
if multifragmentation is at a fixed freeze-out volume, the 
hot system will sample a range of 77/s corresponding to 
different values of pressures at different temperatures. 
While if the breakup volume varies, the different values 
of rj/s found at V/Vq — 3 and 6 in Fig. \5\ implies that 
the system will sample a large range of rj/s in the (P, V) 
plane to give an average. On the other hand, if freeze-out 
is reached at a fixed pressure 23|, |24| (and close to Pc) 
this would result in a rapid increase in ij/s at \T — Tc\ > 
comparable to the rise observed for II2O [14]. 

In summary, we have studied the thermodynamic and 
transport properties in fragmentation of a nucleus within 



a statistical multifragmentation model. For the equili- 
brated system of fragments and light particles evolving 
under Coulomb repulsion, we find the shear viscosity to 
entropy density ratio rj/s exhibits a minimum at a tem- 
perature of T ss 5 MeV in the coexistence phase of in- 
termediate mass fragments and light particles. The mini- 
mum value of (?7/s)min — 2 in multifragmentation is com- 
parable to that at the critical point for liquid-gas phase 
transition in H2O. The temperature dependence of r]/s 
is somewhat sensitive to the surface effects on symme- 
try energy and depends rather strongly on the freeze-out 
volume. 
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